The air shower maximum probed by Cherenkov effects from radio emission 
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Abstract 

Radio detection of cosmic-ray-induced air showers has come to a flight the last decade. Along with the experimental 
efforts, several theoretical models were developed. The main radio-emission mechanisms are established to be the 
geomagnetic emission due to deflection of electrons and positrons in Earth's magnetic field and the charge-excess emission 
due to a net electron excess in the air shower front. It was only recently shown that Cherenkov effects play an important 
role in the radio emission from air showers. In this article we show the importance of these effects to extract quantitatively 
the position of the shower maximum from the radio signal, which is a sensitive measure for the mass of the initial cosmic 
ray. We also show that the relative magnitude of the charge-excess and geomagnetic emission changes considerably at 
small observer distances where Cherenkov effects apply. 
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(— I 1. Introduction 

I Recent experiments (LOPES [2 [5], CODALEMA [3 

O 13]) showed that radio emission can be used as a new com- 
■^-^ plementary technique for measuring cosmic-ray-induced 
j^ air showers. This gave a new impulse to the field of radio 
I— ' emission from air showers and triggered several new experi- 
, ments at the Pierre Auger Observatory [SI [71 [H] , and LO- 
FAR [3] . From the theoretical point of view several models 
have been developed to give a description of the radio emis- 
sion from air showers, such as REAS3 TO;, MGMR [11], 
C^ and more recently ZHAires [T^ , SELFAS 13J , CoReas [H] , 
^ and EVA [15]. 

There used to be a huge discrepancy between the differ- 
ent macroscopic and microscopic models. The microscopic 
models in time domain like REAS [16] were based on the 
• • geosynchrotron emmission of electrons and positrons gy- 
_ J^ rating in Earth's magnetic field. The macroscopic models 
S^ in time domain like MGMR [TB] were based on the macro- 
scopic geomagnetic current due to deflection of charged 
particles in Earth's magnetic field and the net negative 
electron excess in the shower front. Where the micro- 
scopic models in time domain used to predict uni-polar 
pulses, the macroscopic models predicted bi-polar pulses 
with a difference in magnitude up to a factor ten. This has 
only recently been resolved with the addition of a missing 
radiation component to the microscopic models |17j . 

It follows that similar results are obtained for macro- 
scopic models like MGMR [TH] and models using a mi- 



> 



O 
en 






Email address: krijndevries@gmail.com (Krijn D. de Vries) 



croscopic description by adding the electric field contri- 
butions of the separate electrons like REAS3 [lOj. From 
these results it can be concluded that the emission mecha- 
nisms are understood. The main geomangetic mechanism 
is due to the drift of charged particles in Earth's mag- 
netic field. This gives rise to a net macroscopic current in 
the direction of the Lorentz force acting on the particles. 
The time variation of the induced vector potential gives 
rise to a net electromagnetic pulse. A second emission 
mechanism is due to the net negative charge in the shower 
front. The polarization of the signal as a function of the 
observer position with respect to the impact point of the 
shower can be used to disentangle the geomagnetic emis- 
sion and the charge-excess emission [TS] . Furthermore, the 
Lateral Distribution Function (LDF), given by the radio 
signal strength as a function of antenna distance to the 
shower axis, can be used to give information about the 
composition of the cosmic-ray spectrum at the highest en- 
ergies [T8] . 

In [12] it was shown that even though the deviation of 
the index of refraction from unity, the refractivity, is small 
(0(10^'^)), this deviation gives rise to Cherenkov effects 
at realistic observer distances. The importance of these 
effects were discussed in [12l [TH [2Ql |2T| for a simplified 
shower geometry, where the LDF was shown to have a dis- 
tinct peak. In [T3] we showed the first results of the EVA- 
code based on realistic charge and current distributions 
in the air shower which are obtained from Monte-Carlo 
simulations. 

In this article we show that due to Cherenkov effects it 
is possible to accurately extract from the radio signal the 
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position where the shower profile reaches its maximum. 
This is closely related to the mass of the initial cosmic ray. 
In addition Cherenkov effects have a large effect on the rel- 
ative magnitude of the charge-excess emission with respect 
to the geomagnetic emission. The calculations presented 
are done using the EVA |15| simulations to which we give 
a short introduction. 

2. The model 



where we defined the coordinates rf^ = x-^ + r and A = 
hk—h. Here hk is defined as the critical distance behind the 
shower front where Cherenkov effects occur and the vec- 
tor potential diverges. Looking at Eq. Q we see that by 
making use of the finite extent of the particle distributions 
in the shower a non-singular expression is obtained. The 
derivatives of these distributions, w' = dw/dh, and w* = 
dw/dx^, are finite where the vector potential, Ap^, has a 
square-root divergence which is now safely integrated. 



The basis of the EVA model lies at the Lienard-Wiechert 
potentials from classical electrodynamics [32]. The par- 
ticle distributions in the shower front are introduced by 
the weight function 'w{r, h). The total vector potential is 
given by the convolution of the point-like potential and 
the weight function, 



A^^{t,S) = 



d^r 



dhw{r,h)A>'p^{t,x^O, (1) 



where the point-like vector potential given by. 
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is based on the four-current 



4^{t\x:) = jf'{i!)6\x-at'))- 



(2) 
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The vector potential is to be evaluated at the negative re- 
tarded shower time t' . We define S,{t') — —ct'e^w, where 
e^w is the unit vector pointing along the shower axis. The 
velocity is defined asV = c~^d^/dt' . The four- vector R is 
defined by the light-cone condition R^ = L = c{t — t') and 
W = —Ld/d^L. Where L is the optical path length be- 
tween the source ^{t') and the observer position x. In [TH] 
it was already shown that the denominator in the vec- 
tor potential, RV, can be linked to the derivative of the 
retarded time with respect to the observer time RV = 
R^ /{dt' /dt). Furthermore, it was shown that for a real- 
istic index of refraction the denominator could hit a zero 
leading to a divergence in the vector potential. This di- 
vergence is due to the well known Cherenkov effect. To 
overcome these divergences a partial integration is done 
such that the derivatives acting on the denominator \RV\ 
are shifted to the weight function w{r,h). The final ex- 
pression for the electric field is now given by |15| . 
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3. The emission mechanisms 

The main emission mechanism is due to the deflection 
of electrons and positrons in Earth's magnetic field. The 
induced drift gives rise to a net macroscopic current in the 
direction of the Lorentz force. The electric field is given 
by the time derivative of the induced vector potential and 
is polarized in the direction of the Lorentz force. There 
is an additional emission mechanism due to a net electron 
excess in the air shower front. This electron excess induces 
a scalar potential, and the electric field is given by the spa- 
tial derivatives acting on the potential. The polarization 
is given by the spatial derivative and hence becomes ra- 
dial in the direction of the observer with respect to the 
shower axis [23]. Evidence for the charge-excess compo- 
nent in air showers has been observed at the Pierre Auger 
Observatory [21] and at the CODALEMA site [15]. Apart 
from an additional test of the Monte Carlo shower simu- 
lation, the importance of a thorough understanding of the 
charge-excess emission at distances where Cherenkov ef- 
fects play a role lies in the fact that this is the main emis- 
sion mechanism in dense media where Cherenkov effects 
are dominant. This emission is thus crucial for experi- 
ments searching for GZK neutrino's through radio emis- 
sion from showers induced in ice |26j or moon rock [27]. 

The maximum for the geomagnetic radio emission oc- 
curs when the shower maximum is observed at the Cherenkov 
angle. The position of the shower maximum Xmaxis,/ cnj?) , 
is defined by the depth at which the total electromagnetic 
component reaches its maximum. For charge-excess the 
maximum in the radio signal is seen when Cherenkov emis- 
sion is observed from the maximum of the electron excess 
profile. The electron excess profile does not necessarily 
peak at the same height as the total shower profile and 
the Cherenkov peak might be observed at a different ob- 
server distance. Therefore, polarization analysis gives ad- 
ditional information about the relative peak position of 
the total electromagnetic component with respect to the 
charge-excess profile. 

For the charge-excess analysis in [24l [28] , the variable 
R defined as 



R 
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(5) 



is used. By definition the x direction is given by the pro- 
jection of the geomagnetic component eg^o = — e^ x cb 



on the ground plane, where e^ is the unit vector pointing 
along the shower axis and eg is the unit vector pointing 
in the direction of Earth's magnetic field. The y direction 
is given perpendicular to x in the ground plane. With this 
definition R vanishes independent of the observer position 
if there would be no electric field component beside the ge- 
omagnetic emission. Since the charge-excess polarization 
is pointing radially inward from the observer position to 
the shower axis, R vanishes if the observer is positioned on 
the positive x axis, becomes negative moving to the y axis, 
goes to zero again on the —x axis, and becomes positive 
on the —y axis. 

In the following, we consider the simplified geometry 
of a perpendicular incoming air shower with respect to 
Earth's surface. The magnetic field is pointing to the 
North. The observer angle is defined from East (-0 = de- 
grees) to North (-0 = 90 degrees). For the considered ge- 
ometry a sinusoidal pattern should become visible if R is 
plotted as a function of the observer angle ip. In Fig.[T]the 
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Figure 1: The average value of R as defined in Eq. ([5]) 
as a function of tp the observer angle in degrees for a set 
of 24, 10"'^'' eV, showers with a proton as primary shower- 
inducing particle. The primary inducing shower particle 
moves perpendicular to Earth's surface. The squares give 
R for an index of refraction equal to unity, the circles give 
R for a realistic index of refraction following the law of 
Gladstone and Dale given in Eq. (|6|. 

average value of R for a set of 24, 10^^ eV, showers with a 
proton as primary shower-inducing particle is shown. This 
is done for an index of refraction equal to unity as well as 
a more realistic index of refraction following the law of 
Gladstone and Dale, ugd, given by 



nGD = I + 0-226^ pih), 



cm^ 



(6) 



where p{h) is the air density at an atmospheric height h. 
The observer distance d = 75 m is chosen in such a way 
that the shower maximum is observed at the Cherenkov 



angle for geomagnetic emission. For an index of refraction 
of unity we clearly observe a sinusoidal-like pattern. Since 
Cherenkov effects do not affect the polarization of the sig- 
nal we do not expect the qualitative features of this pattern 
to change significantly for a realistic index of refraction. 
This is indeed the case as can be observed from Fig. [T] 
Nevertheless, even though the pattern does not change 
significantly, the maximum value of R increases. 

This is also seen in Fig. [2] where we plot the absolute 
value of R at an observer angle of 90 degrees as a function 
of distance from the shower axis, |_R|(90°; d). The charge- 
excess emission is radially polarized and should thus vanish 
at the shower axis. This is not the case for the geomagnetic 
emission which is polarized in the direction of the Lorentz 
force. It follows that charge-excess emission, and thus R, 
scales with the opening angle with respect to the point of 
maximum emission. For an index of refraction equal to 
unity this is the opening angle with the shower maximum, 
which steadily increases with observer distance, and hence 
\R\ also steadily increases with observer distance. For a 
realistic index of refraction, however, the opening angle 
of maximum emission is equal to the Cherenkov angle and 
stays constant over the range of distances where Cherenkov 
effects are dominant. This is observed in Fig. [2] in the re- 
gion below 140 meters. Below 50 meters, \R\ becomes con- 
siderably larger when Cherenkov effects are included into 
the simulation. At slightly larger observer distances, but 
still in the regime where Cherenkov effects are important, 
\R\ becomes smaller when Cherenkov effects are applied. 
In this region something more subtle is going on. 

This can be understood from Fig. l3J where the shower 
profile (full line) for the total electromagnetic component 
is given as well as the charge-excess profile (dashed line). 
It follows that the charge-excess profile peaks closer to 
Earth's surface (^ 740 g/cm^) than the shower profile (~ 
680 g/cm^). Due to this effect, at 140 meters Cherenkov 
effects have diminished for the charge-excess component, 
but arc still large for the geomagnetic emission and hence 
R becomes small. 

If the charge-excess profile would have peaked at a 
larger height compared to the shower profile for the total 
electromagnetic component, |i?| as plotted in Fig. [2] would 
always be larger when Cherenkov effects are included into 
the simulation compared to the calculations done for an 
index of refraction equal to unity. This is a large effect, 
which should be measurable. The value of \R\ in the regime 
where Cherenkov effects apply is thus sensitive to the rel- 
ative position of the charge-excess maximum with respect 
to the maximum for the shower profile given by the total 
electromagnetic component. 

For distances larger than 140 meters, Cherenkov effects 
become small and the value of R for a realistic index of 
refraction converges to the value of R for a constant index 
of refraction equal to unity. 
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Figure 2: The average value of \R\ as defined in Eq. ^ at 
an observer angle ip = 90 degrees as a function of distance 
from the shower axis d for a set of 24, 10^^ eV, showers 
with a proton as primary shower-inducing particle. The 
shower-inducing particle moves perpendicular to Earth's 
surface. The squares give R for an index of refraction 
equal to unity, the circles give R for a realistic index of 
refraction following the law of Gladstone and Dale given 
in Eq. ([gI 



4. The Lateral Distribution Function: 
The Cherenkov ring 

Due to the extremely thin verticle particle distribution 
close to the shower axis the electric field contains very high 
frequency components when the shower maximum is seen 
at the Cherenkov angle [15]. Away from the Cherenkov 
angle the emitted radiation is not seen at a single instant 
any more and the projection of the shower profile is the 
determining length scale. This is considered as the 'nor- 
mal' varying current radiation scaling with the derivative 
of the shower profile and is limited to lower frequencies. 
This trend for the different frequency components is also 
seen in Fig. |4j Here the Lateral Distribution Function 
(LDF) is shown. The LDF is given by the intensity of the 
radio signal divided by the bandwidth as a function of dis- 
tance to the shower axis. In Fig. [4J the LDF is shown for 
different frequency bandwidths. The simulation is done 
for a typical 10^^ eV proton shower. The LDF at low fre- 
quencies is naturally sensitive to emission on a rater long 
time scale, which is generally given by the 'normal' vary- 
ing current emission and peaks close to the shower axis. 
For the higher frequencies the LDF becomes more sensi- 
tive for Cherenkov emission from close to the shower max- 
imum and peaks further outward. In the 120 — 200 MHz 
band, Cherenkov emission and normal emission compete 
with each other giving rise to a two peak structure. 
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Figure 3: The shower profile given by the full electromag- 
netic component (full line), and the charge-excess profile 
(dashed line) as a function of depth along the shower tra- 
jectory X(g/cm^). The vertical lines denote the position 
of the profile maximum. The charge-excess profile peaks 
closer to Earth's surface than the shower profile given by 
the full electromagnetic component. 
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Figure 4: The Lateral Distribution Function (LDF) for 
different frequency bandwidths. The LDF is given by the 
intensity of the radio signal divided by the bandwidth as 
a function of distance to the shower axis. The LDF is 
calculated for a typical 10^^ eV shower induced by a proton 
primary particle moving perpendicular to Earth's surface. 



5. Determining Xmax from the radio signal 

For an index of refraction equal to unity it is shown 
that one can distinguish between different shower-inducing 
primary particles [TSl |3D]. This is due to the fact that 



at observer distances close to the shower axis the elec- 
tric field is sensitive to the particle distributions in the 
shower front which do not differ significantly for proton 
and iron induced showers. At a large observer distance, 
on the other hand, the electric field is mainly sensitive to 
the shower profile, and thus also to X,„a2:(g/cni^), defined 
by the position where the full electromagnetic component 
of the shower reaches a maximum. By taking ratios at 
different observer distances a handle on Xmax, which is 
closely related to the mass of the initial cosmic ray, is ob- 
tained. 

In this section we will show a similar ratio as well as 
a new, more direct procedure to yield accurate informa- 
tion about X,nax when Cherenkov effects are taken into 
account. 

In [29], it was noticed that the observer position d = 
y/{x-^^Y + {x-^'^Y that corresponds to the Cherenkov an- 
gle can be linked to the emission height along the shower 
axis by. 
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(7) 



for a constant index of refraction n. A more realistic model 
for the index of refraction is given by the law of Gladstone 
and Dale given in Eq. ([6]), where the index of refraction 
depends on the air density and thus atmospheric height. 
To test if the relation in Eq. Q can be used when the 
law of Gladstone and Dale is applied, we simulated a set 
of 100 showers with a 10^^ eV proton as primary shower- 
inducing particle, and a set of 20 showers with iron as 
primary shower-inducing particle. The primary shower- 
inducing particle moves perpendicular to Earth's surface. 
For each of these simulations we determine the peak of the 
LDF, dp, and the position of the shower maximum, Xm, 
which can be easily linked to X^ax- In Fig- [Sl we plot 
Xmax as a function of dp for different frequency band- 
widths. Showers for which the radio signal peaks below 30 
meters have been omitted from Fig. [5] and the following 
analysis. From Fig. [5] it follows that X„iax can be param- 
eterized as 



Xri 



b ■ d„ 



(8) 



A fit with Eq. (IS]) is sufficient to obtain X^ax with an ac- 
curacy of 10 — 15 g/cm^ for the higher frequency bands 
(> 200 MHz). It has been tested that this is independent 
of the energy of the incoming primary particle. In Fig. [6] 
this is shown in the 200-500 MHz band for three different 
energies of 10^^, 10^*, and 10^^ eV. For each energy a set 
of 20 showers with a proton as primary shower-inducing 
particle and a set of 20 showers with an iron atom as pri- 
mary shower-inducing particle is simulated. The fit pa- 
rameters a, and b, and standard deviations are given in 
Table 1. From Fig. [5] it is clear that for the low fre- 
quency bands (< 200 MHz), the sensitivity for X^ax as 
a function of distance to the shower axis is small. For 
the 120 — 200 MHz band a jump is observed in Xmax as 
a function of distance to the shower axis. This jump is 
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Figure 5: The value of Xmax {g/cm^) as a function of dp(m) 
the position of the peak in the LDF for different frequency 
bands for a set of 100 proton (circles) and 20 iron (squares) 
induced showers. 
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Figure 6: The value of A'ma2:(g/cm^) as a function dp(rn) 
fitted by Eq. ^ in the 200-500 MHz band. For three 
different energies, 10^^ eV, 10^^ eV, and 10^^ eV, a set of 
20 proton and 20 iron induced showers has been simulated. 
The standard deviation of the fit is 14.2 g/cm^. 



related to the two peak behavior in Fig. |4] and is due to 
the interplay between the normal radiation and Cherenkov 
emission. Following [181 130] we plot in Fig. [7J Xmax as a 
function of Q300A00 '^hich is defined as 



r,BW 

^d2/dl 



P^^{liV^/m^/MHz; d^d2m) 
P^^iliV^/m^/MHz; d = dlm)' 



(9) 



Table 1: The fit parameters to Eq. 8 for the high frequency bands (> 200 MHz) 
(< 200 MHz). The standard deviation of the fits is also given. 



and Eq. 10 for the low frequency bands 



Band- Width 


scale parameter 


slope parameter 


standard deviation (g/cm^) 


30-80 MHz 


a=161.1 


r;=88.2 


a = 14.4 


120-200 MHz 


a=-208.1 


77=219.9 


cr:=10.2 


200-500 MHz 


a=916.4 


b=4.6 


0- = 14.2 


500-1000 MHz 


a=898.2 


b=3.7 


0- = 11.8 


1000-2000 MHz 


a=860.3 


b=2.9 


G = 12.0 



the power in a fixed Band- With (BW) observed at distance 
d = d2 m from the shower axis divided by the power at a 
distance d = dl m from the shower axis. This analysis is 
limited to bands at low frequencies (< 200 MHz). One of 
the simulated showers peakes to close to Earth's surface 
and has been excluded from the analysis. The simulations 
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Figure 7: The value of Xma2;(g/cm^) as a function Q300/100 
defined by Eq. ^ for the 30-80 MHz band. For three 
different energies, 10^^ eV, 10^^ eV, and 10^^ eV, a set of 
20 proton and 20 iron induced showers has been simulated. 
The standard deviation of the fit is 14.4 g/cm^. 



can be fitted with a logarithmic function given as 



-^Tti.n.T. 



= a-77log(gf2%i) 



(10) 



obtaining an accuracy of ^ 10 — 15 g/cm^. The fit param- 
eters a, and 77, and the standard deviations are given in 
Table 1 for the different low frequency bands (< 200 MHz). 

It should be noted that the obtained results are for a 
simplified geometry of showers moving perpendicular to 
Earth's surface. Also instrumental errors have not been 
included into the analysis. In comparison, the Fluores- 
cence detection technique employed at the Pierre Auger 
Observatory ;3T], determines the shower maximum with 
an approximate accuracy of 20 g/cm^ . 



6. Conclusions 

In this article we have shown that radio emission from 
cosmic-ray-induced air showers can be used to obtain rele- 
vant information about the emission mechanisms and the 
composition of the cosmic-ray spectrum. The influence of 
Cherenkov effects on the emission mechanisms, the geo- 
magnetic emission due to the deflection of charged leptons 
in Earth's magnetic field and the charge-excess emission 
due to a net negative electron excess in the shower front, 
is investigated. This is studied by means of the R param- 
eter which is a measure for the charge-excess component 
in the emission. It is shown that since the polarization is 
not affected by Cherenkov effects, the qualitative behav- 
ior of i? as a function of observer angle does not change. 
Nevertheless, the behavior of i? as a function of observer 
distance to the shower axis is shown to be sensitive to 
Cherenkov effects. It is shown that these effects can be 
used as a direct measure of the relative position of the 
maximum number of particles, X„^ax^ for the full shower 
profile given by the total electromagnetic component with 
respect to the charge-excess profile given by the electron 
excess in the shower front. 

Cherenkov effects determine the structure of the lateral 
distribution function (LDF) defined by the signal strength 
as a function of distance to the shower axis. For the high 
frequency bands (> 200 MHz), it is shown that the peak 
position in the LDF can be linked to X^ax with an ac- 
curacy of 10 — 15 g/cm^. For the low frequency bands 
(< 200 MHz), we use the power ratio, QfJ/^i' given by 
the power at distance d1 divided by the power at distance 
d\ as a measure for Xmax- With this method a similar 
accuracy of 10 — 15 g/cm^ is obtained. These values are 
obtained for a simplified geometry of vertical air showers, 
and instrumental errors have not been included into the 
analysis. 
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